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Lecture  2 

GOVERNING  EQUATIONS,  ASYMPTOTICS,  AND  DEFLAGRATIONS 

The  problem  of  formulating  the  governing  equations  of  combustion  consists, 
at  its  simplest,  in  characterizing  the  flow  of  a  viscous,  heat-conducting 
mixture  of  diffusing,  reacting  gases.  This  is  a  formidable  task  that 
could  fill  a  week  of  lectures  by  itself,  most  of  which  would  not  be  of 
great  interest  to  a  mathematical  audience.  Mindful  of  this,  we  shall  limit 
ourselves  to  a  description,  rather  than  a  derivation,  of  the  simplest 
equations  that  can  be  brought  to  bear  on  combustion  problems.  Only  the 
most  important  assumptions  normally  used  to  Justify  the  equations  will 
be  discussed;  for  a  more  extensive  treatment  the  reader  is  referred  to 
Buckmaster  &  Ludford  (1982}cV»-£) 

We  shall  then  outline  the  asymptotic  method  on  which  the  whole  theory 
rests  and  use  it  to  solve  the  basic  problem  of  combustion:  the  steady 
propagation  of  a  plane  deflagration  wave. 


3p/3t  +  7’(py)  »  0 


(1) 
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I 


where  p  is  the  density  and  y  the  velocity.  But  only  for  dilute  mixtures 
is  the  overall  momentum  balance  identical  to  that  for  a  single  homogeneous 
fluid,  namely 

p  Dy/Dt  *  V-E,  ’  (2) 

where 

E  =  -(p  +  |k7-v)I  +  <[Vy+(7y)T]  (3) 


and  p  is  pressure;  bulk  viscosity  has  been  neglected. 

A  single-fluid  equation  for  energy  balance  is  also  justified,  provided 
account  is  taken  of  the  release  of  heat  by  chemical  reaction.  But  here 
additional  approximations  are  made,  based  on  the  observation  that  tempera¬ 
tures  are  high  and  gas  speeds  low  for  a  large  class  of  combustion 

phenomena  (excluding  detonation);  more  precisely,  a  characteristic  Mach 

-2  -3 

number  is  small  (typically  in  the  range  10  -10  ) .  Then  the  only 

significant  form  of  energy,  other  than  that  of  chemical  bonding,  is  thermal; 
kinetic  energy  makes  a  negligible  contribution  to  the  energy  balance.  For 
the  same  reason,  the  conversion  of  kinetic  energy  into  thermal  energy  by 
way  of  viscous  dissipation  can  be  ignored.  Thus,  when  variations  of  the 
specific  heat  c^  with  temperature  are  neglected,  we  have 


p  CpDT/Dt  -  7* ( XVT)-3p/3t  *  q. 


(U) 


where  q  is  the  heat  released  per  unit  volume  of  the  fluid  by  chemical 
reaction;  the  form  of  q  is  considered  later. 

In  addition,  the  assumption  of  small  Mach  number  implies  that  spatial 
variations  in  pressure  are  small,  so  that  3p/3t  in  the  energy  balance  (U) 
is  due  to  imposed,  uniform  pressure  variations.  We  shall  assume  that  the 
imposed  pressure  is  constant,  i.e.  the  term  vanishes.  The  pressure  term 
in  the  momentum  equation  (2)  cannot  be  neglected,  however;  the  small 
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spatial  variations  are  needed  to  account  for  changes  in  the  weak  velocity- 
field.  A  further  consequence  of  the  virtual  constancy  of  the  pressure 
is  that  the  equation  of  state  of  the  mixture  is  Charles's  law 

pT  =  mpc/R  (5) 

if  the  inert  is  a  perfect  gas.  Here  m  is  the  molecular  mass  of  the  inert, 
Pc  the  imposed  constant  pressure,  and  R  the  gas  constant. 

Consider  now  the  individual  components  of  the  mixture,  denoting  the 
density  of  the  i-component  by  pY^,  where  Y.^  is  the  mass  fraction  and 
i  =  1,2,...,N.  The  reactants  and  their  products  are  convected  with  the 
gas  speed  v,  diffuse  relative  to  the  inert  diluent,  and  are  consumed  or 
generated  by  reaction.  The  diffusion  laws  of  general  mixtures  are 
complicated,  involving  a  diffusion  matrix;  but  for  dilute  mixtures  the 
matrix  is  diagonal  insofar  as  the  reactants  and  products  are  concerned, 
so  that  we  may  write 

pDY./Dt  -  V-Cik.VYJ  =  p.  for  i  =  1,2,...,N-1.  (6) 

Here  p^  is  the  mass  production  rate  per  unit  volume  of  the  ith  component; 
its  precise  form  is  considered  below.  The  equation  for  the  mass  fraction 
YN  of  the  inert  is  more  complicated;  but  it  can  be  obtained  from  the 
relation 

N 

I  Y,  -  1  (T) 

i=l  1 

instead,  once  the  other  Y^'s  have  been  determined. 

Coupling  between  the  fluid-mechanical  equations  (l),  (2)  and  the 
thermal -chemical  equations  (b),(6)  occurs  because  of  density  variations. 

If  these  variations  are  ignored,  the  former  may  be  solved  for  v,  which 
can  then  be  substituted  into  the  latter,  a  substantial  simplication. 
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Such  a  procedure  is  justified  if  the  heat  released  by  the  reaction  is  small, 
but  this  is  not  a  characteristic  of  combustion  systems ,  whose  main  purpose 
is  to  liberate  heat  from  its  chemical  bonds.  For  this  reason,  the  simplified 
system  of  equations  should  be  thought  of  as  a  model  in  the  spirit  of  Oseen's 
approximation  in  hydrodynamics.  However,  to  emphasize  the  mathematically 
rational  nature  of  the  procedure,  we  shall  refer  to  the  simplified  system  as 
the  constant-density  approximation  rather  than  me 'el.  Phenomena  whose 
physical  Dasis  is  truly  fluid-mechanical  (e.g.  the  Darrieus-Landau  instability 
discussed  in  lecture  5)  are  not  encompassed  by  this  approximation,  but 
much  of  importance  is;  it  will  play  a  central  role  in  our  discussion. 

There  remains  the  question  of  the  contribution  of  the  individual 

reactions  to  the  heat  release  q  and  the  production  rates  p  .  It  is 

i 

possible,  in  principle,  to  consider  all  the  reactions  that  are  taking  place 
between  the  constituents  of  a  mixture.  However,  this  is  done  but  infrequently; 
often  a  complete  chemical -kinetic  description  (i.e.  how  the  rates  depend  on 
the  various  concentrations  and  temperature,  or  even  whether  a  particular 
reaction  takes  place)  is  not  available.  Even  when  it  is,  its  complexity 
may  deter  solution  by  anything  short  of  massive  use  of  computers.  For 
these  reasons,  simplified  kinetic  schemes  are  normally  adopted  which  model, 
in  an  overall  fashion,  the  multitude  of  reactions. 

The  simplest  are  the  one-step  irreversible  schemes  that  account  for 
the  consumption  of  the  reactants,  here  taken  to  be  just  a  fuel  and  an 
oxidant.  If  the  reactants  are  simply  lumped  together  as  a  single  entity, 
the  scheme  is  represented  by 

[Y1]  -*■  products,  (8) 

where  brackets  denote  a  molecule  of  the  component  whose  mass  fraction  is 
enclosed.  On  the  other  hand,  if  the  separate  identities  of  the  fuel  and 
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oxidant  are  recognized,  we  have 

v^Y^J  +  \>2[Y2]  Products'»  (9) 

here  the  are  stoichiometric  coefficients,  specifying  the  molecular 

proportions  in  which  the  two  reactants  participate.  We  shall  adopt  the 
scheme  (3)  when  discussing  premixed  combustion  and  (9)  for  diffusion  flames 
(lectures  8-10  only).  These  terms  are  defined  at  the  beginning  of  section  U. 
If  DL  is  the  number  density  of  the  i^*1  component,  so  that 

p .  =  m.N. ,  (10) 

ill 

where  nu  is  the  molecular  mass  of  the  ith  component,  the  reaction  rate  ui  is 
defined  by  the  formula 

N.  =  -v.u.  (11) 

i  l 

It  is  then  common  to  write 

v  61 

id  =  k(T)pT  II  Y  J  (y»0.  positive  consts.  )  (12) 

i  0  J 

for  the  reaction  rate,  an  empirical  formula  that  is  suggested  by  a  theoretical 
treatment  of  so-called  elementary  reactions.  The  product  contains  a  single 
term  for  the  scheme  (8),  two  terms  for  (9)-  The  Arrhenius  law 

k  =  BTae"E,/BT  (B,a,E  consts.),  (13) 

which  we  shall  adopt,  is  at  the  heart  of  our  mathematical  treatment ■,  E  is 
called  the  activation  energy. 

The  heat  release  q  is  a  consequence  of  the  difference  between  the 
heats  of  formation  of  the  products  and  those  of  the  reactnats,  so  that  it 
is  p-oportional  to  u.  Combustion  is  inherently  exothermic,  so  that  we  shall 


write 
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<1  =  Qu.  (lU) 

where  Q(>0)  has  the  dimensions  of  energy. 

2.  Nondimens ional  Equations  and  Their  Shvab-2eldovich  Formulation. 

We  shall  take  units  as  follows: 


temperature  Q/c  E  v  m.  (summation  over  1  or  2  reactants)  (15) 

P  j  J  J 


pressure  p  ,  density  p  ,  mass  flux  M  ,  speed  M  /p  , 
c  r  r  r  r 


(16) 


length  X/c  M  ,  time  Xp  /c  M2,  pressure  variations  M^Vp  .  (17) 

p  r  27  p  r  r  r 

Appropriate  choices  for  the  reference  density  p  and  the  reference  mass  flux 

r 

will  be  made  according  to  the  problem  considered*  The  governing  equations 
in  nondimensional  form  are 


pT  *  mp  c  Z  v  a  /P  RQ,  3p/3t  +  7*(py)  «  0, 
c  P  j  J  J  r 

pDy/Dt  =  -Vp  +  P[v2y  +  |v(V*v)], 

pDT/Dt  -7T  =  JJ,  pDY./Dt  -LT^Y.  =  a.Q, 

l  ill’ 

where  i  runs  from  1  to  N-l,  and 


(18) 

(19) 

(20) 


P  -  kc^/X  (Prandtl  number),  L  =  X/u^c^  (Lewis  number),  (21) 

—A /T  ®  i 

a.  =  -v.m./  Ev.m.  (with  E  a.  =  -l),  0  =  Ve  H  Y,  .  (22) 

1  1 1  j  ■>  J  i  1  j  J 

8  •  Be  Z  v  m  /0.B,  5  •  DM"2,  D  •  [X3q“o'7c1*'‘(I  v  m  )°"1]o,Ta.  (23) 

P  j  J  J  r  p  j  J  J 

Except  in  the  case  y  3  a,  the  so-called  Damkohler  number  D  is  variable; 
it  is  called  a  number  in  spite  of  having  the  dimensions  of  In  the 

context  of  activation-energy  asymptotics  (9  ■*  “),  only  its  value  at  a  fixed 
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temperature  T#  plays  a  role,  so  that  it  may  be  considered  an  assigned 
constant. 

When  one  of  the  Lewis  number  L  is  equal  to  1,  the  differential 
operator  in  its  equation  (20b)  is  identical  to  that  in  equation  (20a). 

We  may  therefore  write 

(pD/Dt  -'72)(T-Yi/a.)  =0,  (2U) 

of  which  one  solution  is 

T  -  Y./a.  5  H.  25) 

li  l 

constant.  If  this  solution  is  appropriate  for  the  problem  at  hand,  Y^  may 
be  eliminated  in  favor  of  T,  thereby  reducing  the  number  of  unknowns.  The 
linear  combusion  (25)  is  known  as  a  Shvab-Zeldovich  variable;  it  is  easier 
to  find  by  virtue  of  satisfying  the  reactionless  equation  (2U). 

3.  Activation-Energy  Asymptotics 

In  these  lectures  we  shall  discuss  a  variety  of  combustion  phenomena 
on  the  basis  of  equations  (l8)-(20).  To  do  this  we  need  an  effective  tool 
for  dealing  with  the  highly  nonlinear  reaction  term  0.  Activation-energy 
asymptotics,  used  in  an  ad  hoc  fashion  by  the  Russian  school  (notably  Frank- 
Kamenetskii  and  Zeldovich)  in  the  kO's,  exploited  in  the  framework  of  modern 
singular  perturbation  theory  (but  in  a  very  narrow  context)  by  aerothermo- 
dynamicists  in  the  60' s,  and  systematically  developed  by  Western  combusion 
scientists  in  the  70 's,  is  just  such  a  tool. 

The  limit  9  -*■  «  is,  by  itself,  of  little  interest:  the  definition 
(22b)  3hows  that  0  vanishes.  To  preserve  the  reaction,  it  is  necessary 
for  V  to  become  unboundedly  large,  i.e.  we  must  consider  a  distinguished 
limit  characterized  essentially  by 
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V  -  e6/T*  (26) 

where  T#  is  a  constant  that  may  have  to  be  found.  The  consequences  of  this 

limit  then  depend  on  the  relative  magnitudes  of  T  and  T#. 

For  T  <  T#,  the  reaction  term  n  vanishes  to  all  orders;  this  is 

B. 

known  as  the  frozen  limit.  For  T  >  T„,  equations  (20)  imply  H  Y  J  -*■  0 

J  J 

exponentially  rapidly,  so  that  at  least  one  Y^  vanishes  like  that,  and 
again  0  vanishes  to  all  orders;  the  equilibrium  limit  hOids.  For  T 
no  more  than  0(9  "*")  from  T*.,  reaction  takes  place,  usually  in  a  thin 
layer  called  a  flame  sheet.  Thus,  with  a  few  exceptions,  the  general  feature 
of  high  activation  energy  is  the  absence  of  chemical  reaction  from  most  of 
the  combusion  field,  the  description  of  which  is  thereby  simplified.  Reaction 
occurs  only  in  thin  layers  (spatial  or  temporal),  whose  description  is  also 
relatively  simple. 

k.  Plane  Deflagration  Waves. 

We  are  now  ready  to  demonstrate  the  efficacy  of  the  technique  that  is 
the  central  theme  of  these  lectures,  by  examining  the  fundamental  problem 
of  premixed  combusion  -  the  plane  unbounded  flame.  But  first  we  pause  for 
a  few  words  on  terminology. 

In  general,  two-reactant  flames  can  be  classified  as  diffusion  or 
premixed.  In  a  premixed  flame  the  reactants  are  mixed  and  burn  when  the 
mixture  is  raised  to  a  sufficiently  high  temperature.  In  a  diffusion  fiame 
the  reactants  are  of  separate  origin;  burning  occurs  only  at  a  diffusion- 
blurred  interface. 

Both  kinds  of  flames  can  be  produced  by  a  3unsen  burner  (figure  1). 

If  the  air  hole  is  only  partly  open,  so  that  a  fuel-rich  mixture  of  gas  and 
air  passes  up  the  burner  tube,  a  thin  conical  sheet  of  flame  stands  at  the 
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mouth;  this  is  the  premixed  flame.  Any  excess  gas  escaping  downstream 
mixes  by  diffusion  with  the  surrounding  atmosphere  and  burns  as  a 
diffusion  flame. 

Separate  origins  do  not  guarantee  a  diffusion  flame,  however.  In 
figure  2  the  reactants  are  originally  separated  by  the  splitter  plate 
but  mix  before  igniting.  The  flame  spread  across  the  oncoming  flow  is 
therefore  premixed.  Behind  this  premixed  flame  the  remaining  portions  of 
the  reactants  are  separate  again,  so  that  a  diffusion  flame  trails  downstream. 

The  plane  unbounded  flame  of  premixed  combustion,  the  so-called  (plane) 
deflagration  wave,  propagates  at  a  well-defined  speed  through  the  fresh 
mixture  and,  accordingly  can  be  brought  to  rest  by  means  of  a  counterflow. 

It  is  natural  to  take  the  mass  flux  of  this  counterflow  as  the  representative 
mass  flux  ;  the  counterflow  is  not  known  a  priori ,  but  is  to  be  determined 
during  the  analysis  of  the  combustion  field.  Indeed,  its  determination  is 
the  main  goal  of  the  analysis.  A  choice  must  also  be  made  of  the  reference 
density  p^;  we  shall  take  it  to  be  that  of  the  fresh  mixture. 

The  continuity  equation  (l3b)  integrates  to  give 

pv  =  1,  ( 27 ) 

so  that 

dT/dx  -  d2T/dx2  *  Q.  (28) 

Since  there  is  only  one  reactant  we  shall  drop  the  subscript  1.  For  1=1 
in  the  corresponding  equation  (20b),  the  Shvab-Zeldovich  formulation  applies, 
showing  that 

T+Y  =  H  *  T^+Y  ST.,  (29) 

ff  o 

where  the  subscript  f  denotes  the  fresh  mixture  at  x  =  (Actually,  H 

is  the  total  enthalpy  of  the  mixture).  Thus, 
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M  *  Z?(Tb-T)e'9/T  (30) 

if  the  most  common  choice  6  =  1  is  made.  Equations  (28)  and  (30)  form 

a  single  equation  for  T,  which  must  satisfy  the  boundary  condition 

T  ■*  as  x  -<=.  (31) 

The  requirement  that  all  the  reactant  be  burnt  provides  the  final  boundary 
condition 

T  -*■  T.  as  x  -*•  +®.  (32) 

b 

Note  that  neither  the  equation  of  state  (l8a)  nor  the  momentum  equation  (19) 
has  been  used;  the  former  provides  p  once  T  has  been  found,  and  the  latter 
then  determines  p  from  v  *  1/p . 

It  is  immediately  apparent,  since  fi  does  not  vanish  for  T  =  Tf, 
that  the  problem  for  T  cannot  have  a  solution.  The  mixture  at  any  finite 
location  will  have  an  infinite  time  to  react  and  so  will  be  completely  burnt. 
This  cold-boundary  difficulty,  as  it  is  known,  is  a  result  of  idealizations 
and  can  be  resolved  in  a  number  of  ways:  the  mixture  can  originate  at  a 
finite  point;  an  appropriate  initial-value  problem  can  be  defined,  without 
the  solution  having  a  steady  limit  of  the  kind  originally  sought;  or  a 
switch-on  temperature  can  be  introduced  below  which  fl  vanishes  identically. 
It  is  one  of  the  virtues  of  activation-energy  asymptotics  that  it  makes 
such  resolutions  unnecessary.  Reaction  at  all  temperatures  below  T* 
(including  Tf)  is  exponentially  small,  so  that  it  takes  an  exponentially 
large  time  for  it  to  have  a  significant  effect;  in  other  words,  T#  is  a 
switch-on  temperature.  Consequently,  it  is  not  necessary  to  discuss  the 
cold -boundary  difficulty  any  further. 

We  now  seek  a  solution  that  is  valid  as  9  •*  *».  Our  construction 
will  be  guided  by  the  assumption  that,  in  the  limit,  reaction  is  confined 
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to  a  thin  sheet  located  at  x  -  0.  On  either  side  of  this  flame  sheet, 
equation  (28)  simplifies  to 

dT/dx  -  d2T/dx2  *  0,  (33) 

which  only  has  a  constant  as  an  acceptable  solution  behind  the  flame  sheet 
(x  >  0),  exponential  growth  being  excluded.  The  boundary  condition 
(32)  then  shows  that  * 

T  =  T.  for  x  >  0;  (3*0 

b 

T,.  is  called  the  adiabatic  flame  temperature.  It  follows  that  the 
b 

temperature  at  the  flame  sheet  is  T,  ,  so  that  this  is  also  the  value 

b 

of  T#  needed  to  specify  the  distinguished  limit  (26) . 

Ahead  of  the  flame  sheet,  equation  (33)  has  the  solution 

T  =  T«  +  Y„ex  for  x  <  0,  (35) 

I  i. 

satisfying  the  boundary  condition  (31)  and  malting  T  continuous  at  x  =  0. 
;;o  structure  could  be  found  for  the  flame  sheet  if  the  temperature  were 
discontinuous. 

Turning  now  to  the  structure ,  which  must  determine  the  still  unknown 
M  (i.e.  V  ) ,  we  note  that  the  form  of  fl  restricts  the  variations  in 
T  to  being  0(e*"1).  Since  the  temperature  gradient  must  be  0(l)  to 
effect  the  transition  between  the  profiles  (3*0  and  (35),  the  appropriate 
layer  variable  is 

5  *  9x}  (36) 

coefficients  in  the  layer  expansion 

T  *  Tb  -  e_1T2d+...  with  <f>  *  {lfT)1  (37 ) 

are  now  considered  to  be  functions  of  equation  (28)  then  shows  that 
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d24>/d52  =  V$e~^  with  5  =  Ve  b/82,  (33) 

while  matching  with  the  solutions  (3*0,  (35)  gives  the  boundary  conditions 

4>  =  -YfC/T2  +  o(l)  as  C  $  *  o(l)  as  £  -*■  +«.  (39) 

In  order  for  this  problem  (in  4»)  to  make  sense,  V  must  be  0(l).  Then 
t>e~9/,~  is  0(82)  in  x  >  0,  so  that  equation  (28)  is  unbalanced  unless  ' 

Y  =  0  (to  all  orders)  behind  the  flame  sheet,  consistent  with  the  result 
(34).  Thus,  equilibrium  prevails  in  x  >  0  even  though  the  mixture  is 
no  hotter  than  the  flame  sheet  there. 

Integrating  equation  (38a)  once,  using  the  condition  (39b),  gives 

( d$ /d§ ) 2  »  29 [ 1- ( 6+1 ) e-^ ] .  (40) 

The  remaining  boundary  condition  will  then  be  satisfied  only  if 

9  «  Y®/2tJ,  (41) 

corresponding  to  the  determination 

_  ?  -9/2T 

M  =  /2D  t5  e  /  Y9  (42) 

r  o  i 

of  the  burning  rate,  (if  D  is  temperature  dependent  it  must  be 
evaluated  at  the  temperature  T^. ) 

Determination  of  the  wave  steed  M  /p  is  the  main  goal  of  the  analysis, 

*  r  r 

and  rightly  so.  But,  at  the  same  time,  the  structure  of  the  combustion 
field  is  obtained  (figure  3).  The  reaction  zone  appears  as  a  discontinuity 
ir.  the  first  derivatives  of  T  and  Y,  a  reflection  of  the  delta  function 
nature  of  9.  in  the  limit  9  +  «>.  Ahead,  the  temperature  rises  and  the 
reactant  concentration  falls  as  the  reaction  zone  is  approached  through 
the  so-called  preneat  zone.  It  is  the  preheat  zone  that  delimits  the 
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combustion  field  and,  therefore,  defines  the  thickness  of  the  flame. 
According  to  the  formula  (35),  sore  than  99?»  of  the  increase  in 
temperature  from  T-  to  T  is  achieved  in  a  distance  5,  i.e.  5X/c  M 

*  D  p 

in  dimensional  terms.  With  this  definition  of  flame  thickness,  we  find 
that  hydrocarbon-air  flames  are  about  0.5  mm  thick.  The  thickness  of 
the  reaction  zone,  which  is  scaled  by  0~\  is  typically  10  or  2C  tines 
smaller. 


5.  Generalizations. 

If  the  reaction-zone  structure  itself  is  required,  equation  (Uo)  must  be  inte¬ 
grated  to  obtain  $  as  a  function  of  £.  The  constant  of  integration  is  fixed 
by  the  boundary  condition  (39a),  which  is  thereby  used  a  second  time. 

The  three  requirements  that  the  boundary  conditions  impose  on  the  solution 


of  the  second-order  equation  (33a)  are  responsible  for  a  definite  value 
of  V.  The  term  laminar-flame  eigenvalue  is  often  used. 

The  analysis  in  section  1*  was  carried  through  with  Levis  number 
unity  for  the  sake  of  simplicity  alone.  ?or  L  ^  1,  the  formula  (1*2) 
is  replaced  by 


M 

"r 


2 

b 


-0/2T 

b/Yfe, 


(U3) 


The  only  change  is  the  replacement  of  *^2D  by  /2LD. 

The  rate  at  which  the  mixture  burns  is  extremely  sensitive  to  the 

_1  o  ,  i 

flame  temperature.  If  7  changes  to  T  -8  T  is,  i.e.  by  an  0(8  ) 

D  0  0 

amount,  the  burning  rate  changes  to 


(1*!*) 


i.e.  by  an  0(1)  amount.  This  result  has  general  validity,  i.e.  it  holds 
whatever  the  nature  of  the  perturbation,  steady  or  unsteady.  The  reason 
is  that,  in  determining  the  eigenvalue  V,  the  perturbation  only  intrudes 
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through  the  matching  of  $  at  £  =  +•  (which  leads  to  the  exponential 
factor);  matching  ahead  of  the  flame  sheet  is  with  an  C(l)  gradient  that 
is  unaffected  by  the  perturbation. 

Questions  involving  deflagration  waves  are  essentially  free-boundary 
problems:  the  location  of  the  flame  sheet  and  the  combustion  fields  on 
either  side  have  to  be  determined  simultaneously.  In  that  context,  the 
structure  problem  (38),  (39),  modified  for  0(0  ^)  perturbations  in 
flame  temperature,  provides  a  leading-order  jump  condition  on  the  normal 
derivative  3T/3n,  namely 

6[(3T/3n)2]  =  Y2  e  **  (U5) 

when  V.  is  taker,  to  have  the  value  (b2).  This,  the  true  role  of  the 
structure  problem,  is  obscured  by  the  analysis  of  the  steady  plane  wave 
giver,  in  the  last  section,  where  M  was  taken  to  be  the  constant  (unknown) 
burning  rate.  If  M  had  been  given  the  value  (12)  without  explanation, 
and  M  used  to  denote  the  (dimensionless)  burning  rate,  then  x  would 
have  been  replaced  by  Mx  in  the  formula  (35)  and  the  jump  condition 
would  have  yielded 

M  =  1,  (U6) 

i.e.  the  location  of  the  flame  sheet. 

Vnen  the  temperature  gradient  vanishes  to  leading  order  behind  the 
flame  sheet,  the  jump  condition  (bp)  gives  the  gradient 

-0./2 

3T/3n  =  Y„e  (U?) 

ahead.  This  result  will  be  needed  repeatedly  later. 

In  many  circumstances,  locating  the  flame  sheet  is  relatively  easy. 

In  the  final  lecture  ve  shall  discuss  several  problems  in  which  it  is  the 


principal  question. 


•V* ' , 
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Figure  Captions 

2.1  Bunsen  burner. 

2.2  Combustion  of  initially  separated  reactants. 

2.3  Profiles  of  T  and  Y  ,  drawn  for  L  =  1  and  T  »  0.25  , 

Yf  =  0.75  .  f 
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